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ABSTRACT 



We study steady and spherically symmetric outflows of pure electron-positron pair plasma as 
a possible acceleration mechanism of relativistic jets up to a bulk Lorentz factor of greater than 
10. We assume that at the inner boundary a "Wien fireball" is realized, which is optically thick 
to Compton scattering but thin to absorption and in a Wien equilibrium state between pairs and 
photons at a relativistic temperature. As was shown by approximate treatments in our previous 
paper, the Wien fireball results in a relativistic outflow by thermal expansion, and thus problems 
with pair annihilation and radiation drag can be avoided. In this paper we present numerical 
solutions obtained with a Monte Carlo simulation of radiative transfer in a relativistic flow. 
Compton scattering, pair annihilation, and pair creation processes are considered in simulating 
the photon trajectories, and we evaluate the photon distribution function, pair creation rate, and 
radiative force. The dynamics of the outflow of pairs are consistently solved with radiative force 
and pair processes by iteration. The numerical results basically confirm our previous finding 
that the formation of powerful relativistic outflows can be obtained by the Wien fireball. Pair 
plasma is relativistically accelerated, and the radiative force does not work as a drag force but 
as an accelerating one because of the relativistic beaming effects. Radiation emitted from the 
photosphere should be observed as MeV peaked emission at infinity with a luminosity on the 
order of the kinetic power of jets. 



Subject headings: elementary particles-galaxies: 
numerical-radiative transfer 



active-galaxies: jets-hydrodynamics-methods: 



1. INTRODUCTION 

The production and bulk acceleration of rela- 
tivistic jets in active galactic nuclei (AGNs) is one 
of the most challenging problems in astrophysics. 
The basic features to be explained are highly rel- 
ativistic velocity with a bulk Lorentz factor above 
10, a huge kinetic power that is almost compara- 
ble to the Eddington luminosity, and collimation 



into a small opening angle (e.g., Ostrowski et al. 
1997; Begelman, Blandford, & Rees 1984). Al- 
though many ideas have been proposed ranging 
from hydrodynamic and magnetic accelerations to 
a radiative one, each model has difficulties for the 
explanation of all aspects of the powerful relativis- 
tic jets. Thus, there is no consensus on how jets 
are produced and accelerated. 

In our previous work (Iwamoto & Takahara 
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2002), we proposed a thermal mechanism of bulk 
acceleration of powerful relativistic outflows on the 
assumption of a spherical symmetric and station- 
ary flow. We assumed that at the base of the 
outflow, electron-positron pair plasma is optically 
thick to Compton scattering but thin to absorp- 
tion. High-energy photons and pairs are assumed 
to be in a Wien equilibrium state. The dynamical 
features of an optically thick pair plasma are quite 
different from those of the optically thin one con- 
sidered in the radiative acceleration model (Phin- 
ney 1982; Inoue & Takahara 1997). We named 
this initial condition of a pair plasma a "Wien fire- 
ball." In the radiative acceleration model, it is well 
known that radiation drag effect (Phinney 1982; 
Sikora et al. 1996) prevents the pairs from obtain- 
ing a bulk Lorentz factor above 5, that the greater 
part of the pairs are annihilated on the way to 
bulk acceleration, and that most of the radiative 
power escapes without scattering. Recently, Be- 
loborodov (1999) studied the behavior of electron- 
positron outflows generated by photon-photon col- 
lisions above the gamma-ray-emitting accretion 
disks. The resultant flow was only mildly rela- 
tivistic because of the radiation drag and Comp- 
ton cooling due to coexisting soft photons in the 
cylindrical geometry. In contrast, our Wien fire- 
ball model supposes that pairs at relativistic tem- 
peratures are generated in a compact region with 
greater compactness and that the effects of the ex- 
ternal soft photons can be neglected in the spher- 
ical geometry. Thus the difficulties of radiation 
drag and pair annihilation are successfully avoided 
owing to the strong beaming of radiation and pair 
production by accompanying high-energy photons. 
It was shown that the pair plasma is relativisti- 
cally accelerated by thermal expansion and that 
most of the pairs survive to infinity. 

However, in our previous paper we made sev- 
eral approximations for a simplified treatment of 
the radiative transfer. There we treated optically 
thick and thin (to scattering) regimes separately. 
We assumed that in the optically thick regime pho- 
tons and pairs form a single fluid while in the opti- 
cally thin regime photons are free-streaming. We 
also assumed that the scattering cross section is 
given by the Thomson one and that the pair cre- 
ation rate is negligible in the optically thin regime. 
In this paper we examine more exactly the features 
of the relativistic outflow that is initiated by the 



Wien fireball, by employing the Monte Carlo treat- 
ment of radiative transfer. In the Monte Carlo 
treatment, a Klein-Nishina cross section is em- 
ployed without using the Thomson approximation 
and the radiation field is faithfully treated both in 
the optically thick and optically thin regimes. In 
addition, pair creation and annihilation are also 
treated by the Monte Carlo method. On the Basis 
of the Monte Carlo simulation, we validly calcu- 
late the angular distribution and energy spectrum 
of photons, the pair creation rate, and the radia- 
tive force. 

In § 2 we describe basic features of the Wien 
fireball model. The difference between the "fire- 
ball" model in gamma-ray bursts (GRBs) and the 
Wien fireball model is noted. In § 3 we present 
basic equations to solve the outflow of the pair 
plasma and the Monte Carlo treatment of the ra- 
diative transfer. In § 4 we show an example of 
the numerical results of pair plasma outflow and 
radiative properties, and in § 5 we present the nu- 
merical results for various boundary values. Fi- 
nally, in § 6 we discuss several issues related to 
the Wien fireball model. 

2. BASIC FEATURES OF WIEN FIRE- 
BALL MODEL 

We assume that a pure electron-positron pair 
plasma, which is optically thick to Compton scat- 
tering, is produced at the base of jet. Such a 
pair plasma will be accompanied by high-energy 
photons because electron-positron pairs are con- 
sidered to be mainly produced by photon-photon 
collisions. Through scattering and pair processes, 
pairs and photons are assumed to be in a certain 
equilibrium state called the "Wien equilibrium" 
(Svensson 1984). 

2.1. Wien Fireball Model for AGN Jets 

For relativistic bulk acceleration, a huge amount 
of energy that exceeds the rest-mass energy is re- 
quired per particle. As such, the fireball model is 
widely accepted for the bulk acceleration model 
of cosmic GRBs (Rees & Meszaros 1992). The 
fireball is a compact high-entropy plasma that 
consists of photons, electron-positron pairs, and 
a small amount of baryons. It is optically thick 
both to scattering and to absorption. It is in a 
complete thermal equilibrium, and the bulk accel- 
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eration should be achieved by thermal expansion. 
In the situation of a fireball, electron-positron 
pairs, photons, and baryons are dynamically cou- 
pled together and there is no need to worry about 
the radiation drag effect because they expand to- 
gether as a single fluid. In the fireball model of 
GRBs, electron-positron pairs are almost wholly 
annihilated in the course of the thermal expansion. 
This is because the temperature decreases below 
the rest-mass energy of the electron before the 
electron-positron pair creation/annihilation pro- 
cesses become frozen. When the thermal expan- 
sion lasts until all the internal energy is expended, 
the attainable bulk Lorentz factor of the outflow 
is simply given by r m <~ E/(m p N), which repre- 
sents the mean free energy per baryon normalized 
by its rest mass m p . Here E and N denote the to- 
tal energy flux and total number flux of baryons, 
respectively, and we use the unit of c = l(c is 
the speed of light). For steady, spherically sym- 
metric, and relativistic flows, the bulk Lorentz 
factor and temperature behave as Y = To(r/ro) 
and 9 = 9o(r/r )~ 1 , as considered in the fireball 
model of GRBs (Goodman 1986; Paczynski 1986). 
Here r is the radius of the fireball, 9 denotes the 
temperature normalized by electron mass, and the 
subscript denotes the physical quantities at the 
inner boundary. 

Can this fireball model be directly applied to 
AGN jets? There is no possibility of a fireball 
forming because the characteristic size of AGNs is 
much larger than that of GRBs. The characteristic 
size of AGNs is 3r g ~ 10 14 cm, where r g = 2GM 
is the gravitational radius for the black hole mass 
M (the typical value of M is taken to be 10 8 M ), 
with G being the gravitational constant. If a com- 
plete thermal equilibrium is assumed, the temper- 
ature of the radiation is estimated as about 10 5 
K, which is far below the rest-mass energy of the 
electron, and no pair creation is expected. The 
rest-mass density of baryons in an optically thick 
plasma is far greater than the energy density of 
photons, and no relativistic flow is expected. How- 
ever, if the radiation is optically thin to absorp- 
tion, a complete thermal equilibrium is not real- 
ized and an alternative possibility exists. From 
the size and luminosity of AGNs, a pair plasma 
is expected to be optically thin to absorption but 
thick to scattering. It is possible to form copious 
electron-positron pairs when high-energy photons 



are provided sufficiently. Even for mildly relativis- 
tic temperatures, copious pairs can be produced 
for a finite optical thickness to scattering (see the 
next subsection). In such a state, photons and 
electron-positron pairs are not in a complete ther- 
mal equilibrium, but they are coupled by electron 
scattering and pair processes. In this state, the 
photon density is far below that in the complete 
equilibrium state at a given temperature. The 
physical state of such a pair plasma may be charac- 
terized by a Wien equilibrium state in the simplest 
case. The Wien equilibrium state is known to be 
strongly constrained by thermal balance and pair 
equilibrium owing to an inefficiency of pair an- 
nihilation at relativistic temperatures (Lightman 
1982; Svensson 1984; Kusunose & Takahara 1985). 
Such a Wien equilibrium state may be established 
at the base of AGN jets, which we call a "Wien 
fireball" (Iwamoto & Takahara 2002). Although 
the Wien fireball is not in a complete thermal equi- 
librium, the dynamical feature of the Wien fireball 
is the same as that of the normal fireball in GRBs. 
Pairs and photons are dynamically coupled to each 
other and relativistically expand together as a sin- 
gle fluid. While the pair plasma is optically thick 
at the base, it becomes optically thin after the 
expansion and electron-positron pairs are dynam- 
ically decoupled from photons outside the photo- 
sphere, where pairs and photons do not behave 
as a single fluid any longer. As discussed in our 
previous paper (Iwamoto & Takahara 2002) and in 
subsection 2.3 below, pairs can survive without ex- 
cessive annihilation in this model. In Figure 1 we 
illustrate a schematic diagram of the pair plasma 
outflow. 

2.2. Optical Thickness at the Base of the Jet 

In this subsection we argue that the base of the 
relativistic jet is expected to be optically thick to 
scattering. Let us consider the optical thickness of 
the jet in the compact region where its formation, 
bulk acceleration, and collimation should occur. 
We can easily estimate the optical thickness under 
some basic assumptions. The jet is assumed to be 
composed of almost pure electron-positron pairs 
with a solid opening angle £lj , a bulk flow velocity 
/3, and a Lorentz factor T. The longitudinal optical 
thickness is estimated by using the Thomson cross 
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Fig. 1. — Schematic diagram of a pair plasma flow- 
ing out from the optically thick regime to the thin 
one. In the optically thick regime, pairs and pho- 
tons behave as a single fluid and expand together. 
In the optically thin regime, photons stream out 
almost freely from the photosphere, and pairs and 
photons do not form a single fluid any longer. 
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Here m p and m c denote the rest mass of a pro- 
ton and an electron, respectively, (7c } denotes the 
mean Lorentz factor of random motion of electron- 
positron pairs, n c denotes the total number den- 
sity of electrons and positrons, and r denotes the 
distance from the central black hole. The kinetic 
luminosity of the jet and the Eddington luminosity 
are denoted by Lj and LEdd, respectively. Because 
the term m p /m e is 1836 and because the collima- 
tion factor 47r/f2j is much larger than unity, the 
optical thickness r exceeds 1 when the other terms 
are not much smaller than 1, i.e., when Lj is com- 
parable to ^Edd and when (7 C ) is not excessively 
high for r of within several times r g . Note that T 
is not large at the base of the jet. Therefore, AGN 
jets should be optically thick to scattering at the 
scale of the gravitational radius, where jet forma- 
tion must occur. By setting r = 1 in equation (1), 
the radius of photosphere is also estimated as 
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The characteristic radius of the photosphere r p h 
for AGNs is about 10 - 100r g when T of around 
10 is attained there. This consideration seems to 
be fully consistent with the Wien fireball model, 



in which the bulk Lorentz factor increases in pro- 
portion to r in the optically thick regime. 

2.3. Survival of Electron-Positron Pairs 

For the pair plasma outflow, there has been 
an argument that pairs disappear in the course of 
the outflow (Celotti & Fabian 1993; Blandford & 
Levinson 1995). If the copious electron-positron 
pairs are confined in a compact region, catas- 
trophic pair annihilation should occur. Here we 
briefly discuss why pairs can survive in the Wien 
fireball model in spite of this. Conventionally, 
whether pairs annihilate or survive in the course of 
the expansion may be assessed by comparing the 
annihilation time scale t ann = n e +/(n e +) ann with 
the dynamical one idyn = r /(T/3). Here n e + and 
(n c +) ann are the number density of positrons and 
its annihilation rate, respectively. When t ann is 
shorter than idyn, the pair annihilation is consid- 
ered to be effective, and this is written as 
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which is almost the same as the optically thick 
condition for Compton scattering (n c -o~^R/Y > 
1) unless 5 > 1. Here n c - is the number den- 
sity of electrons, and we use the annihilation rate 
given by Svensson (1982) in his equations(57) and 
(59). This coincidence results from the fact that 
the cross section of pair annihilation and that of 
Compton scattering are the same order. There- 
fore, the condition t ann < idyn is satisfied in the 
optically thick regime inside the photosphere. In 
this regime, sufficient pair creation by high-energy 
photons compensates for pair annihilation as long 
as the temperature remains relativistic (6 > 1). 
Since the time scale of pair creation is comparable 
to that of pair annihilation, the number of pairs 
is not determined by the time scale of pair anni- 
hilation. Not only photons injected at the inner 
boundary but also those produced by pair annihi- 
lation contribute to the creation of pairs, so that 
a quasi-pair equilibrium state is realized. 

On the other hand, when i ann > td yn , the pair 
annihilation is ineffective. The number of electron- 
positron pairs is almost conserved because the pair 
processes are frozen. There is no need to keep 
the temperature relativistic in the optically thin 
regime, so there is no worry about the disappear- 
ance of pairs in the Wien fireball model. In our 
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previous paper, we gave the freezeout radius rf r , 



Under these boundary conditions at r 



ro 



which represents the location at t a 



yn- 



It 



was shown that r& is almost the same as the pho- 
tospheric radius r p h and that the survival of pairs 
is achieved in the case of the "relativistic freeze- 
out" (Iwamoto & Takahara 2002). 

3. FORMULATION 

We assume that pure electron-positron pairs 
and high-energy photons are injected at the in- 
ner boundary r = r . The optical thickness of 
such a region to electron scattering is large, but 
that to absorption is smaller than unity. We as- 
sume for simplicity that pairs and photons are in a 
Wien equilibrium state at a temperature 9q on the 
boundary. We fix the bulk Lorentz factor at the 
boundary as T = and use the total lumi- 

nosity as a free parameter. The optical thickness 
To is calculated from E and 6$. The number den- 
sities of pairs and photons in a Wien equilibrium 
are related by (Svensson 1984) 
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Here Kj is the 7-th order modified Bessel function 
of the second kind. For pure electron-positron pair 
plasmas (n c + — n e -), the number density ratio is 
given by n c +/n 7 = A(9). The total luminosity E 
is written in terms of the number injection rates 
of positrons N e +^, electrons N c - , and photons 
i\T 7i o at the inner boundary (which are integrated 
over the spherical surface 47rr 2 ) as 



E = 



30 o iV 7 ,o + (7o)(A r c+,o + ^e-,o) 



(5) 



Here (70) is the mean Lorentz factor of ran- 
dom motion of pairs at the boundary. Using 
equation(4), we get N c +_ , N e - <0 , and j\T 7i o as 
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■$6 {) + 2( 1{) )A{6 Q )m c ' 
These number injection rates are also represented 
by the number densities n c +,o an d %,o at the 
boundary as 

N c + fi = 47rr 2 n c +, r /3o, = 47rro 2 n 7 , r /3o- 

(8) 



2r g , we solve the dynamics of pairs and radia- 
tion. For the dynamics of pair plasma, we solve 
the equations for a relativistic perfect fluid inter- 
acting with radiation, which are described in sub- 
section 3.1. Electrons and positrons are assumed 
to obey the Maxwell-Boltzmann distribution. For 
the radiative transfer, we use a Monte Carlo sim- 
ulation in which the effects of the Klein-Nishina 
cross section and the anisotropic distribution of 
photons are considered. The method of the Monte 
Carlo treatment is described in subsection 3.2. In 
this treatment we can assess the radiative force 
and the pair creation rate by photon-photon col- 
lision precisely. The effects of nonuniformity and 
the relativistic motion of pairs are also dealt with 
adequately. The radiative force is obtained by 
evaluating the energy and momentum losses of 
photons per unit volume of a fluid element, and 
the pair creation and annihilation rates are given 
by the loss and gain of the number of photons. In 
order to solve the dynamics of pair plasma out- 
flow consistently with the radiative transfer, we 
iteratively solve the dynamics and the radiative 
transfer. At first, we simulate trajectories of pho- 
tons by the Monte Carlo technique for a given flow 
of pair plasma, which is taken from the approxi- 
mate numerical solutions of Iwamoto & Takahara 
(2002). Then we obtain the radiative force and the 
net creation rate of pairs and solve the dynamics 
of pair plasma outflow again. We make several it- 
erations until the solution converges and we finally 
obtain a consistent solution. 

3.1. Dynamical Equations for Pairs 

We treat the electron-positron pairs as a rela- 
tivistic fluid suffering from interactions with pho- 
tons. We consider the energy and momentum 
transfers by Compton scattering and by pair pro- 
cesses. They are represented as "radiative force" 
F^. For a steady and spherically symmetric flow, 
conservation laws of energy, momentum, and num- 
ber of pairs are given by 
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respectively. Here p c and P e are the energy density 
and pressure of electron-positron pairs, which are 
given by 

p c = 2m c n c+ (7), (12) 



P c = 2m c n»+0, 



(13) 



respectively. The total number density of elec- 
trons and positrons is n c + + n e - = 2n c + (since 
n c + = n e - for a pure electron-positron pair 
plasma), and n c + is the net creation rate of pairs, 
i.e., pair creation rate minus annihilation rate. 
The pair annihilation rate (?i c +) alln is calculated 
as (Svensson 1982) 



-cr T n 0+ n e - 



1 



20 2 



ln(1.120+ 1.3) 



l 

(14) 



where v denotes the thermal velocity of pairs. The 
right-hand side of the energy and momentum con- 
servation laws denotes the radiative force in the co- 
ordinate frame, which is given through the Lorcntz 
transformation of the radiative force in the comov- 
ing frame of a fluid element. The pair creation rate 
and radiative force are calculated by the Monte 
Carlo simulation described in the next subsection. 

Although we solve the optically thick and thin 
regimes with a single method, the notion of the 
photospheric radius r p h, where the optical thick- 
ness to the scattering becomes unity, may be im- 
portant to understand the numerical results. It 
should be noted that r p h depends on the energy 
and direction of the photon and that most pho- 
tons are moving outward in the radial direction 
in the relativistic outflow. Although the energy 
dependence of the cross section is significant for 
photons around MeV range, we define the pho- 
tospheric radius r p h conventionally by using the 
Thomson cross section as 

/>oo 

t= dr (n e + +n e -)<r T r(l -(3) = 1. (15) 

Because of the Klein-Nishina effect, the real pho- 
tosphere is located at a somewhat smaller distance 
than this conventional one. 

3.2. Monte Carlo Treatment of Radiative Trans- 
fer 

The Monte Carlo method for Compton scat- 
tering was described in Pozdnyakov et al. (1977, 



1983), where the Klein-Nishina cross section is 
considered. Here we should also take into ac- 
count pair annihilation and creation, and we need 
to calculate the distribution of photons, too. In 
addition, we should take the bulk flow into ac- 
count. Because the flow velocity is relativistic 
and nonuniform, we divide the calculation regime 
of pair plasma flow into many spherical shell ele- 
ments (see Figure 2), in practice 3323 shells from 
r = 2r g to 10 4 r g with an equal logarithmic in- 
terval. We set the physical quantities of the pair 
plasma (the bulk velocity T/3, temperature 9, and 
number density n c + [= n e -]) to each fluid element. 

We simulate the trajectory of each photon 
propagating through the pair plasma and pho- 
ton field, where the photon interacts with elec- 
trons, positrons, and photons themselves. Ele- 
mentary processes to be considered arc Compton 
scattering and electron-positron pair annihilation 
and creation. While the pair annihilation pro- 
cess generates photons, the pair creation process 
makes photons vanish, and Compton scattering 
exchanges photons. We simulate all these pro- 
cesses by the Monte Carlo method. In the simu- 
lation, we employ a photon splitting method by 
assigning a weight to a simulated photon and trac- 
ing many trajectories for each injected photon. 

Photons are injected both from the inner 
boundary and from pair annihilation throughout 
the flow. At the inner boundary, we generate pho- 
tons obeying an isotropic Wien distribution in the 
comoving frame (ro = y/3/2) with the same tem- 
perature as the electron-positron pairs. Then, we 
transform them to the coordinate frame. For the 
photons arising from pair annihilation, we gen- 
erate two photons per event. Pair annihilation 
photons are generated according to the rate given 
in equation (14). We randomly choose the anni- 
hilating pairs (electrons and positrons) from the 
Maxwell-Boltzmann distribution by taking into 
account the dependence of annihilation proba- 
bility on the cross section and relative velocity 
between the annihilating particles. We calculate 
the energy and momentum of photons from those 
of annihilated pairs by using conservation laws of 
energy and momentum. The spectrum of the an- 
nihilation photons and the method of calculation 
are taken from Ramaty & Meszaros (1981). 

We generate 01b simulated photons at the inner 
boundary; for most cases we take 9Tb = 2 x 10 5 . 
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Conforming to the number injection rate of real 
photons at the inner boundary N lt o, the weight 
parameter of the simulated photons Wb is assigned 
as 



(16) 



Each photon is given an equal weight at injection, 
for simplicity. 

For the pair annihilation photons we generate 
and inject 29t a = 200 photons per each fluid ele- 
ment (the factor 2 arises from the fact that each 
annihilation event generates two photons) so that 
the total number of simulated annihilation pho- 
tons is 6.646 x 10 5 for 3323 fluid elements. We set 
the weight parameter of the annihilation photons 



w a as 



1 

9tl 



(n e +)ann SV '— = — ()l c +)ann SV. 

(17) 

Here SV and SV are the comoving volume and 
the coordinate volume of the fluid clement, respec- 
tively In equation (17), the factor St' /St converts 
the annihilation rate in the comoving frame to that 
in the coordinate frame. [We note that (?i. +)ann 
means the pair annihilation rate not in the co- 
ordinate frame but in the comoving frame.] At 
the second equality, we use the Lorentz invariance 
St'SV — StSV. These weight parameters are used 
when the physical quantities of the radiation field 
are calculated. The weight parameter is suitably 
divided into those of split-simulated photons so as 
to conserve their sum. The weight parameter de- 
creases after the photon-photon absorption. 

At first, we calculate the probability of interac- 
tion within the shell Pint and that of escaping into 
the adjacent shell P osc . These are given by 



P csc = exp(-Ttot), Pnt = 1 - P ef 



(18) 



where r to t is the total optical thickness to all 
the interaction processes (Compton scattering and 



pair creation; r to t = t c , 



-T C re) along the trajec- 



tory in the shell. Here T CO mp and t cic represent 
the optical thickness for Compton scattering and 
that for electron-positron pair creation to the shell 
boundary, respectively. The simulated photon is 
divided into two parts: the "escaping" part and 
the "interaction" part (see Figure 2). The photon 
weight of each part is divided in proportion to the 
probability P osc and p nt . The trajectory of the 



escaping part is regarded as entering the adjacent 
shell and that of the interaction part suffers from 
further interaction before escape and continues to 
be simulated to trace the trajectories of photon 
again and again. 

For Compton scattering, we use the Klein- 
Nishina cross section <7kn> and T comp is calculated 
considering the Maxwellian distribution function 
/ e of scattering particles (Pozdnyakov et al. 1983) 
as 

vp), (19) 



Al J d 3 pf e a KN (l 



where v is the velocity of the particle, /i is the 
cosine between the directions of the photon and 
particle, and Al denotes the path length of the 
photon trajectory to the boundary of the adja- 
cent shell. For pair creation, the cross section for 
photon-photon collision taken from the rep- 

resentation in Gould & Schreder (1967). Optical 
thickness for pair creation r cro is calculated from 
the distribution function of photons / 7 as 

T crc = Al J dflg J de / 7 (e, fig) cr 77 (e, i, Q)v R . 

(20) 

Here e and e are the energies of target and incident 
photons, fig is the propagation direction of the 
target photon, is the collision angle between the 
target and incident photons, and i>r(/x) = 1— cos 
is the relative velocity between photons. 

Thus, we get the probabilities for Compton 
scattering and pair creation (which means photon 
vanishing) : 



p — 

1 comp — 



7"comp ~l~ 7"cre 



Tcomp ~l~ 7"cre 



"Pnt- 

; (21) 

We choose Compton scattering or pair creation by 
this ratio of probabilities by random number gen- 
eration. When we choose pair creation, we make 
the photon vanish and quit tracking the trajec- 
tory of the photon for this part. When we choose 
Compton scattering, we simulate the Compton 
scattering and further trace the scattered photon 
as in Pozdnyakov et al. (1983). Then we continue 
tracking the trajectory to the next interaction and 
escape. 

At each interaction point, the weight parame- 
ter of the interaction part decreases by a factor 
of p nt . If the weight parameter of the interac- 
tion part becomes e w times smaller than that at 
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Fig. 2. — The schematic picture of the trajec- 
tory of a simulated photon propagating in the pair 
plasma outflow. Since the flow and photon distri- 
bution are highly nonuniform, the spatial region 
divides into many spherical shell elements. 



the injected point, we finish tracing the photon 
trajectory. For the nondimensional technical pa- 
rameter e w (e w < 1), we set at 10~ 12 . By using 
the weight parameter like this, photons that have 
far smaller probability than 1 /91 a or 1 /91b can be 
adequately treated (Pozdnyakov et al. 1977). 

After simulating N photons (typically we take 
N = 2 x 10 5 ), photon distribution in the phase 
space is calculated as a superposition of simulated 
photons. The energy density and flux of photons 
are obtained. We can also obtain radiative force 
and pair creation rate by accumulating the energy, 
momentum, and number losses of simulated pho- 
tons. While the pair annihilation rate is derived 
from the distribution of electron-positron pairs, 
the pair creation rate is derived from that of pho- 
tons as follows. 

(n e+ ) cre = - J dQ, e J dQ s J de / 7 (e) / de / 7 (e)cr 77 (e, e, 

where tt e is the direction of the incident photon. 

4. AN EXAMPLE OF NUMERICAL RE- 
SULTS 

We have performed numerical simulations for a 
wide range of parameters of the total luminosity 
E and the boundary temperature 9q. In this sec- 
tion we present a typical example of the numerical 
results for inner boundary conditions 

E/L Edd = 24.9 , O = 2 , with r = 85.4. (23) 

We first present the results of the pair plasma out- 
flow and next those of the radiation field. The 



number, energy, and momentum transfers from 
the radiation field to the pair plasma are also pre- 
sented. 

4.1. Pair Plasma Outflow 

Here we show the result of the pair plasma out- 
flow. The top panel of Figure 3 shows the behav- 
ior of the bulk velocity Yj3 and the temperature 
9. It is seen that the photosphere is located at 
r p h = 9.98r g and that the temperature rapidly de- 
creases with radius, becoming 6* p h = 0.311 at the 
photosphere. The bulk Lorentz factor at the pho- 
tosphere turns out to be r p h/3 p h = 5.52. These nu- 
merical values agree well with the analytic predic- 
tion of r ph ~ 8.84r g , 6 ph ~ 0.452, and r ph ~ 5.41 
(Iwamoto & Takahara 2002). The behavior in- 
side the photosphere represents the basic feature of 
relativistic thermal expansion, as expected. Out- 
side the photosphere, the pair plasma is further 
accelerated by the thermal energy at the photo- 
sphere and by the additional acceleration by ra- 
diative force. The acceleration of bulk motion is 
saturated at about r ~ 100r ff and the terminal 
Lorentz factor is T = 12.2. Compared with nu- 
merical results in our previous paper, the behavior 
of the outflow is quite similar except for small dif- 
ferences in numerical values. The location of the 
photosphere slightly shifts outward. This leads to 
a slightly smaller 6> p h and a slightly larger r p h- 
There are also some differences in the asymptotic 
quantities. Compared to the analytic estimates, 
the asymptotic bulk Lorentz factor r M is smaller 
by about 25%, the number fraction of electron- 
positron pairs is larger by about 10%, and the ki- 
netic' luminosity of pairs is smaller by about 10%. 

The middle panel of Figure 3 shows the rel- 
ative number fraction of pairs and photons as a 
function of the radial coordinate. The pair frac- 
tion decreases steadily according to the decrease of 
temperature. Because pair creation and annihila- 
tion processes are unfrozen inside the photosphere, 
the decrease of temperature leads to the decrease 
of pair fraction. Outside the photosphere, the de- 
crease of pair fraction is modest and becomes neg- 
ligible for r > 50r g . About 30% of the initial pairs 
annihilate and about 70% survive in this calcula- 
tion. Therefore, this result confirms our prediction 
that the pair annihilation problem can be avoided 
successfully. The bottom panel of Figure 3 shows 
the luminosities of pairs and radiation as a func- 



8 



: \ 






: \ 












10 100 


1000 


100 



1000 10000 




100 1000 10000 

r/r. 



Fig. 3. — Top: Behavior of velocity and tempera- 
ture for E/L E dd = 24.9, O = 2, with t = 85.4. 
The solid line denotes the bulk Lorentz factor of 
the pair plasma and the dash-dotted line de- 
notes the equilibrium Lorentz factor of radiation 
field r cq /3 oq . The dashed line denotes the tem- 
perature of the pair plasma. Middle: Fraction of 
the particle and photon numbers. The solid and 
dashed lines denote the number fractions of pairs 
and photons, respectively. Bottom: Luminosities 
of various components. The dashed line denotes 
the total luminosity of pairs and radiation. The 
solid line denotes the kinetic luminosity of pairs, 
while the dotted line denotes the luminosity car- 
ried in a form of the rest mass of pairs. Thus, the 
interval between the solid and dotted lines denotes 
the luminosity carried in the form of the thermal 
energy of pairs, and that between the dashed and 
solid lines denotes the luminosity carried by radi- 
ation. 



tion of the radial coordinate. The total luminosity 
turns out to remain constant by necessity. Inside 
the photosphere, the kinetic luminosity of pairs 
slightly decreases because pairs are converted to 
photons. Outside the photosphere, the kinetic lu- 
minosity of pairs is kept almost constant, where 
the thermal energy of pairs at the photosphere 
converts to bulk kinetic energy by thermal expan- 
sion. Although pairs lose their kinetic luminosity 
by their annihilation process to some degree, ac- 
celeration by radiative force compensates for the 
annihilation loss of the kinetic luminosity. The 
terminal kinetic luminosity of pairs turns out to 
be 52% of the total luminosity, while the remain- 
ing 48% is carried away by photons. The spectrum 
and the angular distribution of these photons are 
shown in the next subsection. 

4.2. Radiation Field 

4-2.1. Spectrum of Photons 

In Figure 4 we present the spectrum of photons 
at the inner boundary (ro = 2r g ) and at infin- 



ity (r a 



10 T g ), both measured in the coordi- 



nate (observer's) frame. The former is the "in- 
put spectrum," which represents the spectrum of 
injected photons at the inner boundary. The lat- 
ter is the "output spectrum," which represents the 
spectrum of escaping photons at the outer bound- 
ary. Compton scattering and pair annihilation and 
creation change the input spectrum into the out- 
put one. Basically, the output spectrum is not so 
different from the injected one, and the injected 
Wien spectrum of photons is almost conserved. 
This fact can be comprehended by the energy con- 
servation law per photon. Photons (contained in 
the Wien fireball) are at a relativistic tempera- 
ture at the inner boundary, but the bulk motion 
of the photons is only mildly relativistic. As the 
Wien fireball expands, the temperature of photons 
decreases, but at the same time the Lorentz fac- 
tor of the bulk motion increases, which boosts the 
energy of photons by the relativistic beaming ef- 
fect. The effect of decreasing temperature, i.e., the 
mean energy of photons in the comoving frame, is 
compensated for by the relativistic beaming effect. 
In the course of the expansion, the internal energy 
of the Wien fireball is converted to the kinetic en- 
ergy of the bulk motion according to the energy 
conservation law. Thus, the mean energy of pho- 
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Fig. 4. — Spectrum of simulated photons. The 
shaded histogram denotes the spectrum of injected 
photons at inner boundary r — ro = 2r g , which we 
call the "input spectrum." The open histogram 
denotes that of escaping photons from the outer 
boundary r = 10 4 r 9 , which we call the "output 
spectrum." Both spectra are measured in the co- 
ordinate (observer's) frame. Here, e/mc 2 is the 
photon energy normalized by the electron rest- 
mass energy. 



tons in the coordinate frame is almost conserved. 

Two features of the output spectrum are noted. 
First, the number flux of output photons is some- 
what larger than that of the input ones because 
net pair annihilation increases the photon num- 
ber. Second, the output spectrum is somewhat 
broader than the input Wien spectrum. The lat- 
ter is mainly due to Compton scattering in the op- 
tically thin regime, which tends to make the pairs 
and photons depart from an equilibrium. Not only 
photons injected at the inner boundary but also 
those injected by pair annihilation are scattered 
by electrons and positrons. Down-scattering by 
cooled pairs reduces the energy of photons and 
forms the low-energy tail of the output spectrum. 
On the other hand, the high-energy part of out- 
put spectrum also enlarges a little, which is also 
due to Compton scattering. The temperature of 
electron-positron pairs is not uniform, and scatter- 
ing by the pairs at multiple temperatures broad- 
ens the photon spectrum. The energy dependence 
of the cross section (Klein-Nishina regime of the 
cross section) affects the broadening too (at higher 
energies we see a deeper and higher temperature 
region). The output spectrum is peaked around a 
few MeV, and its luminosity is comparable to the 
kinetic luminosity of pairs. 

4-2.2. Angular Distribution of Photons 

In Figure 5 we present the angular distribu- 
tion of photons in the comoving frame of the pair 
plasma. Angular distribution is almost isotropic 
in the optically thick regime. The histograms of 
photons (Figs. 5 a and 56) show an almost flat 
distribution. On the other hand, those in the op- 
tically thin regime show a skew distribution. As is 
seen in Figure 5 c? and 5e, outgoing photons (jj, > 0) 
dominate. Obviously, the degree of skewness in- 
creases with r (Fig. 5a — > bb — > 5c — > 5c? — > be). In 
the optically thin regime (Figs. 5c! and 5e), there 
is a small fraction of ingoing photons that are scat- 
tered back or emitted inward in the optically thin 
part. At the photosphere (Fig. 5c), the angular 
distribution of photons is mildly anisotropic. 

4-2.3. Pair Creation and Annihilation Rates 

In Figure 6 we present the pair creation rate, 
which is calculated from the Monte Carlo simula- 
tion, and the pair annihilation rate, which is cal- 



10 



(a) r = 2r g (b)r = 5r g 




(c) r = 9.98 r g 

0.3 1 — i — i — i — i — i — i — i — I 

0.25 - 
0.2 - 

0.15 - | 
0.1 - J 
0.05 ^^^^B 

-1 -0.5 0.5 I 
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Fig. 5. — Angular distribution of photons at var- 
ious shells (measured in the comoving frame of 
the pair plasma). These histograms represent the 
fraction of the photon number oriented to cos -1 // 
from the radial direction. The radial coordinates 
are (a) r = 2r ff , (b) r = 5r g , (c) r — 9.98r g , (d) 
r = 15r 9 , and (e) r = 100r 9 : (a) and (b) are in 
the optically thick regime, (d) and (e) are in the 
optically thin regime, and (c) corresponds to the 
photosphere. 



culated from the physical quantities of electron- 
positron pairs. The pair annihilation rate de- 
creases with increasing r since the number den- 
sities of electrons n e - and positrons n e + decrease. 
The pair creation rate also decreases with r be- 
cause of the decrease of the number density of pho- 
tons n 7 and the increase of its anisotropy. Since 
we assumed that electrons and positrons obey an 
isotropic Maxwell-Boltzmann distribution, the an- 
nihilation rate is simply obtained by equation (14). 
Let us compare the pair creation rate with the 
pair annihilation one. Inside the photosphere 
(r < r ph), the pair creation rate is almost equal to 
the annihilation one, indicating that pair creation 
balances pair annihilation. On the other hand, 
outside the photosphere (r > r p h) pair creation 
does not balance pair annihilation anymore. The 
pair creation rate is less than a thousandth of the 
pair annihilation rate and is negligible. The main 
reason for this is that the angular distribution of 
photons becomes anisotropic in this regime. 




r/r g 

Fig. 6. — Pair creation and annihilation rates. 
The solid line denotes the pair annihilation rate, 
the triangles denotes the pair creation rate, and r 
is the radial coordinate. 

4-2.4- Radiative Force 

Figure 7 shows the energy transfer rates be- 
tween photons and pairs per unit comoving vol- 
ume. The energy transfers through Compton scat- 
tering, pair annihilation, and pair creation are 
shown separately. The energy transfers for these 
three processes decrease with r because the num- 
ber densities of electrons, positrons, and photons 
decrease. The energy loss rate of pairs by pair an- 
nihilation is almost in proportion to the pair anni- 
hilation rate, and energy gain of pairs by pair cre- 
ation is roughly proportional to the pair creation 
rate. In the optically thick regime, the energy loss 
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and gain almost balance. The energy gain of pairs 
by Compton scattering is smaller than that by pair 
processes in the optically thick regime. In detail, 
the photons slightly give energy to pairs by these 
three processes (F° > 0), but the radiative force 
is not as effective as the gas pressure force of the 
pairs themselves. The pairs and photons are bulk- 
accelerated by their own internal energy, and their 
luminosities do not vary much. In the optically 
thin regime, the energy gain by pair creation be- 
comes smallest because of the steep decrease of the 
pair creation rate. 




Fig. 7. — Energy transfer rates between photons 
and electron-positron pairs. This graph represents 
the energy transfer rates for three processes in the 
comoving frame of the pair plasma. The solid line 
denotes the energy loss of pairs due to pair annihi- 
lation, the open triangles denote the energy gain 
of pairs due to pair creation, and the open circles 
denote the energy gain of pairs through Compton 
scattering. 




r/r 8 

Fig. 8. — Momentum transfer rates from photons 
to electron-positron pairs in the comoving frame 
of pair plasma. The symbols are the same as in 
Figure 7. The open triangles denote the momen- 
tum gain of pairs for pair creation, and the circles 
denote that for Compton scattering. There is no 
momentum transfer for pair annihilation. 



Figure 8 represents momentum transfers through 
pair creation and Compton scattering. The be- 
havior is almost the same as the energy transfers 
shown in Figure 7. Note that there is no momen- 
tum transfer for pair annihilation in the comoving 
frame of a fluid element because electrons and 
positrons are isotropic in the comoving frame. In 
both the optically thick and optically thin regimes, 
electron-positron pairs gain momentum from pho- 
tons by pair creation and by Compton scattering 
(the same as the energy transfer). In the optically 
thick regime, the momentum gain by pair creation 
and that by Compton scattering are comparable 
to each other, but in the optically thin regime, the 
momentum gain by pair creation decreases faster 
because of a steep decrease in the pair creation 
rate. Therefore, the momentum gain by Compton 
scattering dominates and affects the pair plasma 
dynamics. As seen in Figure 8, photons give mo- 
mentum to pairs and accelerate the pair plasma 
further, which indicates that radiation does not 
act as the drag force but does as the accelerating 
one. 

5. NUMERICAL RESULTS FOR VARI- 
OUS BOUNDARY VALUES 

We have also computed for various boundary 
values, and the results are tabulated in Table 1. 
As is seen, the terminal Lorentz factor proves to 
be 5 ~ 20 and the terminal kinetic luminosity of 
pairs accounts for about 40%-60% of the total lu- 
minosity for 1 < 8 < 4 and 0.3 < E/L Edd < 30. 
Here we examine the dependence of the physi- 
cal quantities at the photosphere on To and show 
that they agree well with the analytic predictions 
(r p h oc T 1/3 r , ph oc to~ 1/3 6» , and r ph oc r 1/3 ; 
see Iwamoto & Takahara (2002)). Then we de- 
scribe the properties of the asymptotic values of 
the bulk Lorentz factor and the kinetic luminosity 
of the pair outflows. 

Figures 9, 10, and 11 show r p h, p h, and F p h 
as functions of to, respectively. We see that the 
numerical results agree well with the analytic pre- 
dictions. As is shown in Figure 9, r p h is quite 
tightly correlated with To. In a few cases that 
correspond to the "nonrelativistic freezeout," the 
numerical result deviates from the analytic pre- 
diction because a modest disappearance of pair 
plasma near the photosphere occurs at subrela- 
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tivistic temperatures. Figure 10 shows 8 p h/0 as 
a function of to- Numerical solutions roughly fol- 
low the analytic prediction, but they seem to give 
a slightly lower temperature than the analytic pre- 
diction for larger to. When To is further increased 
for constant E, 6 and 9 p h decreases, leading to 
the nonrelativistic freezeout. A few symbols that 
deviate from the line upward correspond to this 
case, where a lower temperature causes excess an- 
nihilation of pairs inside the photosphere. This 
leads to a smaller r p h and larger 9 p h/0o because 
the spatial interval for thermal expansion becomes 
shorter. Figure 11 indicates that T p h is also tightly 
correlated with To, and numerical results can be 
interpreted in the same way as in Figures 9 and 
10. These features of the photosphere confirm the 
predictions of the thermal expansion inside the 
photosphere of Wien fireball even when we take 
account of detailed processes of radiation transfer 
such as the Klein-Nishina cross section and effects 
of anisotropic photon distribution. 




Fig. 10. — Relation between the photospheric tem- 
perature p h and the initial optical thickness To- 
The solid line represents the analytic prediction 
Oph/^o = t -1 / 3 . The configuration of the sym- 
bols is the same as in Fig. 9. 
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Fig. 9. — Relation between the photospheric 
radius r p h and the initial optical thickness To. 
The solid line represents the analytic prediction 
r P h/ro — to 1 / 3 . The crosses, plus signs, filled 
triangles, open squares, and filled diamonds re- 
fer to the cases E/L^dd — 0.3, 1, 3, 10, and 
30, respectively. The open diamond refers to 
E/L Ed d = 24.9. 

Next, we present the properties of the termi- 
nal outflow. In Figure 12, we present the ter- 
minal Lorentz factor. The terminal Lorentz fac- 
tor mainly depends on the temperature at the in- 
ner boundary, since the electron-positron pairs are 
thermally accelerated both inside and outside the 
photosphere. The terminal Lorentz factor is ana- 
lytically estimated as about 4#o for the relativistic 
perfect fluid. It predicts that a higher boundary 
temperature results in a larger terminal Lorentz 




Fig. 11. — Relation between the photospheric 
Lorentz factor r p h and the initial optical thick- 
ness To- The solid line represents the analytic pre- 
diction r p h/ro = to 1 / 3 . The configuration of the 
symbols is the same as in Fig. 9. 
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factor. In addition, beamed radiation from the 
photosphere further accelerates the pair plasma 
in the optically thin regime. Thus, the terminal 
Lorentz factor is larger than the simple analytic 
prediction of 4#o- For nonrclativistic freezeout 
(# p h < 0.5), this additional acceleration effect by 
radiation becomes large enough and does not 
decrease below 5 even at a nonrclativistic temper- 
ature of #o = 0.5. 

In Figure 13 we present the terminal kinetic lu- 
minosity divided by the total luminosity as a func- 
tion of # p h- This ratio indicates the efficiency of 
energy conversion into the production of matter 
outflow. An efficiency of 2/3 (dotted line) corre- 
sponds to the energy fraction of electron-positron 
pairs in a Wien equilibrium state at a relativistic 
temperature. For relativistic freezeout, the effi- 
ciency turns out to be about 50%-60%, which is a 
little smaller than the canonical value of 2/3. This 
is due to the energy loss for pair annihilation pro- 
cesses. For nonrelativistic freezeout, the efficiency 
becomes small because of the further decrease in 
pairs. The remaining part of the luminosity is ac- 
counted for by the MeV peaked radiation. 
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Fig. 12. — Terminal Lorentz factor Too/^oo. The 
abscissa is the temperature at the boundary #o- 
The configuration of the symbols is the same as in 
Fig. 9. 

Summarizing these results, the attainable 
Lorentz factor of electron-positron pair plasma 
is more than 10 for the inner boundary condition 
of 9q > 2, and the kinetic luminosity is compa- 
rable to the total luminosity, provided that Wien 
equilibrium states of pure electron-positron pairs 
are prepared at the inner boundary with a tem- 
perature of a few times the electron mass and an 
optical thickness to the scattering of more than 5. 
When the photosphcric temperature becomes non- 
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Fig. 13. — Ratio of the terminal kinetic luminosity 
of pairs to the total luminosity. The abscissa is 
the photospheric temperature. The configuration 
of the symbols is the same as in Fig. 9, and the 
dashed line denotes the canonical value of 2/3. 

relativistic because of low boundary temperature 
and/or large optical thickness at the boundary, the 
terminal kinetic luminosity becomes low compared 
with the total luminosity, although the terminal 
Lorentz factor is still relativistic. The differences 
between the present results based on the detailed 
Monte Carlo simulation and the previous results 
based on simplified treatments of radiative trans- 
fer are surprisingly small. Generally, they are less 
than 10%, and at most about 20%-30%. 

6. CONCLUSION AND DISCUSSION 

We have considered relativistic outflow of elec- 
tron positron pairs generated by the Wien fire- 
ball and obtained consistent numerical solutions of 
the dynamics of pairs and radiative transfer. We 
have utilized Monte Carlo simulation for a faithful 
treatment of radiative transfer including Compton 
scattering obeying the Klein-Nishina cross section, 
pair creation and annihilation, and we have ob- 
tained consistent solutions by several iterations. 
Compared with the results in our previous pa- 
per (Iwamoto & Takahara 2002), in which simpli- 
fied radiative transfer was adopted, the differences 
in numerical results are very small, generally less 
than 10%, and at most about 20%-30%. Thus, we 
have confirmed basic features of the Wien fireball 
model and the analytic predictions described in 
our previous paper. In addition, we have obtained 
a detailed energy and angular spectrum of pho- 
tons, which enables us to calculate accurate radia- 
tive force. Basically, our present results indicate 
that our previous method overestimate the radia- 
tive force by about 10%; thus, the correct terminal 



14 



values of the Lorentz factor and kinetic luminosity 
become smaller while that of the number fraction 
of pairs becomes larger, compared with the results 
of the simplified treatments. 

We briefly reiterate the basic features of the 
Wien fireball model. This model successfully 
achieves relativistic bulk acceleration and high ki- 
netic power. It can avoid the problems with pair 
annihilation and radiation drag that have long 
been an annoyance with pair models. At the same 
time, it predicts a strong MeV peaked emission 
with a power similar to the kinetic power. These 
features are mainly due to the assumptions that a 
Wien fireball of electron positron pairs at relativis- 
tic temperature is produced and that the effects 
of external soft photons are neglected. Although 
whether such situations are realized or not is a crit- 
ical issue for the present model, these issues have 
not been quantitatively discussed (see the discus- 
sion in Iwamoto & Takahara (2002)). Our view 
is that instantaneous heating produces a runaway 
production of relativistic pairs due to slow anni- 
hilation and resultant dynamical expansion. The 
radiation field in the fireball is dominated by the 
internally produced gamma rays rather than ex- 
ternal softer photons. For reference, the gamma- 
ray compactness of the fireball Z 7 ~ ^ Tttdd ~Jl ' s 
very large, nominally about 10 4 -10 5 , and about 
10 2 -10 3 even when the collimation is taken into 
account (see below), as supposed. As a result, 
photons and pairs will quickly settle in a Wien 
equilibrium. When the fireball has a conical ge- 
ometry in reality, the external soft photons will 
be reflected on the side of the fireball. The shield- 
ing process discussed by Illarionov & Krolik (1996) 
can help to realize our situation since gamma rays 
from the fireball and X-rays from the outer disk 
make a pair atmosphere on the surface of the fire- 
ball and the atmosphere shields the penetration of 
soft photons into the fireball. 

Since we have assumed spherical symmetry, 
we cannot solve the collimation problem, which 
should be investigated in the future with two- or 
three-dimensional treatments. Within spherical 
symmetry, the total luminosity of the successful 
Wien fireball must be higher than the Eddington 
luminosity. However, the required total luminos- 
ity becomes smaller if the outflow is collimated 
within a small opening angle, as is actually ob- 
served. Even when the jet is collimated to the solid 



angle Qj , the dynamics will not be much different 
from the spherical one and the total luminosity of 
jet becomes smaller by a factor VL-JAt:. Since f2j 
is typically as small as 10~ 2 , the required total 
luminosity can be smaller than the Eddington lu- 
minosity. As was discussed in our previous paper, 
there remain several unsolved issues with the Wien 
fireball model. These include the realization of the 
initial states in the hot accretion disk and its hot 
corona, or other ways to produce copious electron- 
positron pairs, the problem of baryon contamina- 
tion, and collimation into the jet. These issues will 
be investigated in the future, along with observa- 
tional confrontation of the Wien fireball model. 
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E/L Edd 9 t 


r p h r ph /3p h p h 




0.3 0.5 3.23 
0.3 1 1.78 


2.74 1.45 0.367 
2.25 1.23 0.809 


5.15 0.130 0.433 
7.05 0.164 0.547 


1.0 0.5 10.1 
1.0 1 6.63 
1.0 2 3.27 


3.95 2.12 0.264 
3.59 2.18 0.481 
2.68 1.71 1.25 


6.00 0.378 0.378 
6.91 0.525 0.525 
11.3 0.601 0.601 


3.0 0.5 28.4 
3.0 1 19.1 
3.0 2 10.2 
3.0 3 6.53 


5.15 2.66 0.232 
5.20 3.02 0.343 
4.19 2.81 0.765 
3.33 2.34 1.46 


7.63 0.904 0.301 
8.15 1.45 0.483 
10.4 1.74 0.580 
15.3 1.77 0.590 


10 0.5 86.8 
10 1 61.7 
10 2 33.6 
10 3 22.2 
10 4 15.8 


6.88 3.40 0.200 
7.19 3.93 0.285 
7.08 4.26 0.435 
5.76 3.85 0.810 
4.75 3.46 1.30 


9.20 2.27 0.227 
10.6 3.95 0.395 
10.8 5.66 0.566 
13.6 5.62 0.562 
17.5 5.48 0.548 


24.9 2 85.4 


9.71 5.46 0.312 


12.2 13.0 0.522 


30 1 182 
30 2 98.8 
30 3 65.5 


9.91 5.11 0.215 
9.87 5.80 0.311 
8.83 5.64 0.512 


12.1 9.65 0.322 
13.1 15.9 0.530 
14.6 16.5 0.550 



Table 1: Numerical results for various boundary 
conditions. First, second and third columns refer 
to the total luminosity, boundary temperature and 
boundary optical thickness. The fourth through 
sixth columns show the numerical results for the 
quantities at the photosphere, while the seventh 
through ninth columns show the numerical results 
for those at infinity. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 
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